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ABSTRACT 

Context. During the evolution of rotating first stars, which initially consisted of only hydrogen and helium, CNO ele- 
ments may emerge to their surface. These stars may therefore have winds that are driven only by CNO elements. 
Aims. We study weak wind effects (Gayley-Owocki heating and multicomponent effects) in stellar winds of first gener- 
ation stars driven purely by CNO elements. 

Methods. We apply our NLTE multicomponent models and hydrodynamical simulations. 

Results. The multicomponent effects (frictional heating and decoupling) are important particularly for low metallicity 
winds, but they influence mass loss rate only if they cause decoupling for velocities lower than the escape velocity. The 
multicomponent effects also modify the feedback from first stars. As a result of the decoupling of radiatively accelerated 
metals from hydrogen and helium, the first low-energy cosmic ray particles are generated. We study the interaction of 
these particles with the interstellar medium concluding that these particles easily penetrate the interstellar medium of 
a given minihalo. We discuss the charging of the first stars by means of their winds. 

Conclusions. Gayley-Owocki heating, frictional heating, and the decoupling of wind components occur in the winds of 
evolved low- metallicity stars and the solar metallicity main-sequence stars. 

Key words, stars: winds, outflows - stars: mass-loss - stars: early-type - hydrodynamics 
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Many aspects of chemical evolution in our Universe remain 
unclear. The first elements, helium and a very small amount 
of, e.g., lithium, were almost certainly s ynthesised durin g 
the era of primordial nucleosynthesis fe.g.. lCoc et al.|[2004l) . 
The first stars in the Universe can theref ore be considered 
to be purely hydrogen-helium stars (see iLoeb et "aTl l2008l 
for a review). 

Subsequent chemical evolution is less clear, partly be- 
cause it is difScult to test the theoretical predictions obser- 
vationally. For example, stars with a very low abundance 
of iron are observed, which are expected to be relics from 
ancient times (e.g.. IChristlieb et al, 2002i ). However, their 
relevance as observational testbeds to the theory of evo- 
lution of chemical com position may be hamp ered by sev- 
eral secondary effects (jTornatore et al.l^2007^ . Several hy- 
potheses have been developed to explain t he chemical com- 
position of these stellar relics (cf.. Shigevama et aLll2003l 
' lUmeda &: Nomotoll2005llMevnet et al.ll2006[ l. 



Hot star winds are supposed to play an important role 
in the chemical evolution of our Universe. Since they re- 
move material from the outer stellar envelopes, they only 
affect the stellar mass during the early phases of stel- 
lar evolution and do not contribute to the change in the 
chemical composition of the interstellar medium. On the 
other hand, as soon as freshly synthesised elements emerge 
at th e stellar surface during later pha ses of stellar evolu- 
tion (jMevnet et al.ll2006llHirschill2007D . the hot star winds 
may contribute to the chemical evolution of the interstellar 
medium even before the star possibly explodes as a super- 
nova. 



Hot star winds are studied mainly by assuming a so- 
lar mixture of elements and infor mation about winds of 
more exotic composition is scarce ([Vink fc de Koteil I2005L 
iGrafener fc HamannI I2008L IUnglaubll2008^ . An interesting 
mixture of heavier elements, which is uncommon in con- 
temporary Universe, is represented by a pure CNO com- 
position. This composition may be typical of later phases 
in the evolution of the first stars. This is connected to the 
possibility that the envelopes of the first stars in later evolu- 
tionary phases arc enriched by the products of helium burn- 
ing (Mcynct et al. 2006, Hirschi 2007). A chemical mixture 
rich in CNO elements and underabundant in iron is typi- 
cal for one group of low-me tallicity stars (e.g.. iNorris et al.l 
IT997I . IChristlieb et al.ll200^ . 

The study of CNO driven winds is importan t not only 
for early stellar generations (cf.. lUnglaubl l2008l) . The low 
density winds of present stars are also accelerated mostly 
by CNO lines because the cont ribution of other heavier 
elements is relatively small fe.g.. lVink et al.ll2001l) . 

To understand the role of CNO driven winds in 
hot evo lved first stars, iKrticka Sz Kuba3 (|2009L hereafter 
iPaper J ) calculated wind models of these stars. They con- 
cluded that CNO elements do not drive winds as efficiently 
as iron peak elements because of the lower number of their 
strong lines. Therefore, the total amount of mass lost by 
these winds does not significantly affect stellar evolution. 
For subsequent stellar generations, the wind enrichment of 
primordial halos by heavier elements does not overcome 
the metallicity threshold for the formation of very massive 
stars. On the other hand, the enrichment could be large 
enough to change the behavior of primordial stars during 
their formation ~ a mass fraction of CNO higher than about 
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10^^° is sufficient to enable hydrogen burning via CNO cy- 
cle and preclude initial helium burning by means of the 3a 
reaction ( Cassisi & Castellani f993). 

Some CNO driven winds (especially the low-metallicity 
ones) may be subject to weak wind effects. For weak winds, 
two effects that are negligible for high-density winds may 
become importa nt, namely the Gayley-O wocki (Doppler) 



Decome importa nt, namely tne Ijayiey-U ' 
heating/cooling (jGavlev fc Owockil 1994 



hereafter 



and multicomponent effects ( Krticka &: Kubatll200lL here 
after mil. 

The GO heating/cooling is caused by a frequency dif- 
ference between photons entering and escaping the Sobolev 
resonance zone. Multicomponent effects are connected to 
the momentum transfer between heavier elements (acceler- 
ated by line absorption) and bulk wind material, i.e., hy- 
drogen and helium. In low-density winds, momentum trans- 
fer may become inefficient, causing fr ictional heating or 



even decoupling of wind c ompo nents (IC astor et al 



Springmann fc PauldrachI 1992L [k rticka fc Kubat 



Owocki fc Puldl200aiVotruba et al.ll2007llUnglaubll2008D . 



1976 . 



2001 



To understand the role of weak wind effects in the CNO 
winds of massive first stars, we calculate models of the mul- 
ticomponent winds of these stars for which GO heating is 
taken into account. 



2. Model assumptions 

The parameters of the model stars studied (see Table [J) 
were obtained from the grid of the evolution ary models of 
initial ly zero-metallicity stars calculated by iMarigo et al.l 
(|2001f) . The parameters of these model stars were selected 
to cover a large area of the HR diagram. Only stars that 
may have wind for the mass frac tion of h eavier elements 
Z < 10~^ are considered here (see iPaper J ) . 

For the study of a multicomponent stellar wind, we 
used two different codes. We firs t used a NLTE niulti- 
compoii ent stationary wind code (jKrticka fc Kubatll200ll 
iPaper IT ), briefiy described in Sect. 12.11 This code enables 
us to predict a wind mass-loss rate and hydrodynamical 
structure, but does not allow us to study the evolution of 
the decoupled fiow . Second, we used a t ime-dependent mul- 
ticomponent code (jVotruba et al.l[2007t ) to study the decou- 
pling of wind components (see Sect. 12 2p . 



2.1. NLTE multicomponent wind models 

A reader interested in a more detailed description of the 
latest v ersion of our NLTE multicomponent wind code can 
consult IPaper I| for more detailed information on NLTE 
equations and KKIj for description of multicomponent hy- 
drodynamic equations. Here we only summarize the basic 
features of our models. 

To calculate NLTE multicomponent wind models, we 
assume a spherically symmetric stationary stellar wind. The 
excitation and ionization state of considered elements is 
derived from the statistical equilibrium (NLTE) equations. 
The ionic models are either taken from the OSTAR2002 
grid of model stellar atmospheres (jLanz fc Hubenvl [20031 
|2007() or prepared by ourselves. The ionic mod els used here 
are based mainly on t he Opaci ty Projec t data (Peach et al.' 



1988, 



1992 . 



Luo fc Pradhiin 1989. Tullv et al.. 1990. Seaton et aL 
Fernlev et aOigggt ). For more details of our ionic list. 



we refer to iPaoer 



Table 1. Radius i?,, mass M, and effective temperature 
Toff of studied model stars 





Evolved stars 






ZAMS stars 




Model 




M 


Tc« 


Model 


k. 


M 


Tcff 




[R©] 


[Mq] 


[kK] 




[Rq] 


[Mo] 


[kK] 


M999-1 


56.4 


100 


29.9 


M200 


1.65 


20 


65.3 


M999-2 


125 


100 


20.1 


M120 


1.42 


12 


49.9 


M999-3 


510 


100 


10.0 


M070 


1.30 


7 


34.8 


M500-1 


11.1 


50 


50.0 


M040 


1.17 


4 


23.6 


M500-2 


33.7 


50 


29.9 










M500-3 


72.0 


50 


20.6 










M500-4 


303 


50 


10.1 










M200-1 


4.1 


20 


50.0 










M200-2 


19.9 


20 


24.5 










MlOO-1 


11.4 


10 


20.2 










MlOO-2 


45.6 


10 


9.8 










M050-1 


5.1 


5 


20.1 











The solution of the radiative transfer equation is sim- 
plified for both continua (neglect of li ne transitions ) and 
lines (using the Sobolev approximation. [CastorifTOT^ ). The 
line radiative force is calculated in the Sobolev approx- 
imation using NLTE level populations. The Sobolev ap- 
proximation and in particular the critical point approach 
for the c alculation of mass-loss rates ha s gained some 
criticism (iLucvl 120071 iMiiller fc Vi"^l2008ll . However, we 
note that the comoving-frame calculation of the radiative 
force agrees with the Sobolev approximat ion in the super- 
sonic part of smooth line-driven winds ([Pauldrach et al.l 
I1986L iPulsl [19871. Moreover, the mass-loss rates derived 
from hydrodynamical simulations agre e with those derived 
using CAK critical p oint approach ([Owocki et al.l I1988L 
iFeldmeier et al.l Il997f) . Atomic data for the line radia- 
tive force calculations are take n from the VALD database 
Kupka e t al. 19 99), a nd also partly 



1995,, 



(Piskunov et al. ^^ ^-^ ^ ^ 

from IWiese et all ()1996[ ) and [kuruczl (fi^93) . The surface 



emergent flux (i.e., the lower boundary condition for the ra- 
diative transfer in the wind) is taken from the H-He spheri- 
cally symmetric NLTE model stellar atmospheres of lKubatI 
( 200^ and references therein). 

The derived radiative force (including the force due to 
the light scattering on free electrons) is used to solve the 
hydrodynamic equations. We solve the equation of conti- 
nuity, momentum, and energy equations for each compo- 
nent of the flow. To calculate the radiative cooling/heating 
term, we use the thermal balance of electrons method 
([Kubat et al.l Il999f) taking all considered bound-bound, 
bound- free, and free- free transitions into account. 

For our present purposes we calculated five-component 
wind models with wind components corresponding to car- 
bon, nitrogen, oxygen, and free electrons, and a passive 
wind component (hydrogen and helium). The inclusion of 
GO heating into our models is described in Sect. 13.11 

2.2. Time-dependent multicomponent wind models 

To calculate time-dependent models, we restrict ourselves 
to a ID spherically symmetric, isothermal, quasineutral, 
two-component outflow consisting of metals (namely oxy- 
gen, carbon, and nitrogen) that scatter stellar photons 
in numerous spectral lines, and passive plasma (con- 
sisting of hydrogen and helium). We use the simplified 
two-component model instead of a more accurate five- 
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component one because we wish to suppress numerical in- 
stabilities and keep the problem finitely computable. The 
acting forces in our model are gravity, dynamical friction, 
gas pressure gradients, and the line radiative force, the last 
of which acts only on line-scattering ions. 

The continuity equations are (see Eqs. (1) and (2) of 
IVotruba et alll2007[ ) 



1 d{ 



dpp ^ 

dt dr 

dt dr 



and the equations of motion are 



dvp 
dt 



dvi 
'dt 



dvp 

dvi 
' dr 



1 dpp 

Pp dr 

Pi dr 



Pp 



Soft , 



- 5rad - 9cS 



■pi 



(la) 
(lb) 

(2a) 
(2b) 



In these equations, Pi, Wi, and pi represent the density, ve- 
locity, and pressure of metals, respectively, and Pp,Vp, and 
Pp denote the same quantities for passive plasma. The line 
radiative force is denoted by gl^^. The effective gravita- 
tional acceleration is gcs = GM{l — Te)/r'^, where Fg is the 
Eddington factor, i.e., the ratio of radiative force caused by 
electron scattering to gravitational force and G is the grav- 
itational constant. In this two-component model, radiative 
force caused by Thomson scattering acts on both compo- 
nents , absorbing ions and passive plasma (jOwocki &: Puis! 
|2002|) . The frictional force Rpi between metals and passive 
plasma is expressed as 



^^pi — ripriikpiGi^Xpi) . 



(3) 



are the number densities of passive plasma and 
respectively, fcpi is the frictional coefficient. 



where i 

absorbing ions, 

G(a;pi) is the Chandrasekhar function, and X n\ is the dimen- 
sionle ss drift speed between components (see lVotruba et aLl 
[20071 Eqs. (11) and (12)). 



To solve the hydrodynamic equations Eqs^ 



CIS. (HI) and Q 
1 IVotruba et al.l 



we use the hydrodynamic code described in IVotruba et 
( 2007(1 , which provides technical details of the code and de- 
scribes the numerical schemes used. Here we changed the 
method of calculation of the frictional term. Dynamical fric- 
tion is the most difficult term to compute. Due to big stiff- 
ness of the system of partial differential e quations, we ap ply 
a fully implicit scheme for this term (see IVotruball20 1ff). 

We adopt a relatively small Courant number 0.05, which 
can stabilize the numerical problems caused by strong de- 
coupling instability. This instability is caused by the depen- 
dence of the frictional force on the drift velocity between 
spe cies (for more details about the dec oupling instability, 
see lOwocki fc Pulsll2002L iKrticka fc Kubat 20021 As the fi- 
nal Courant time step, we use the minimum of the separate 
time steps calculated for both individual flow components. 



3. GO heating 

3.1. Inclusion in the NLTE models 

The Gayley-Owocki (GO, Doppler) heating term (per 
unit of volume) is given in the Sobolev approxima- 
tion by summation of heating contributions over all 



lines (jCavlev fc Owockil 11994 iKrticka fc Kubatll200ll see 

Eq. (28) therein) 



GO 



2ttv ^ 

2 Z^'"tl,.iVijIcl[tij,(J,p*) 

TC 

lines 



(4) 



where v is the radial wind velocity, the thermal speed of 

1/2 

ion with mass mj is Vth,\ = {'2,kT /m{) , Ic — 4:Hc is the 
core intensity (frequency-dependent), — Xijcr/ {vijv), 



t I dfi 



da; x(j){x) exp 



t^{x) 



the frequency-integrated line opacity is 

Xij 



-] dih,^ 



(5) 



(6) 



where n^, Uj, gi, and gj are number densities and statisti- 
cal weights of levels giving rise to the line with oscillator 
strength fij and frequency i/y , D{iJ,) is unity for p, > p^ and 

1 /2 

zero otherwise (//* = (l — R-t/T"^) )j core penetration and 
escape probabilities are given by 



1 1 — exp 



t 



and 



/3 



1 1 — exp 



( L- 



dp 



dp, 



(7a) 



(7b) 



respectively, (/)(a;) is the line profile (assumed to be given 
by a Gaussian function). 



^{x) = / dx'(l){x'), (8) 

J X 

and the variable a was introduced bv ICastorl (|l974f ) to be 



r dv ^ 
V dr 



(9) 



This variable determines the sign of the GO heating. For 
cr > 0, the sign of the function j(t, a, /i*) is negative, caus- 
ing the wind GO cooling (typically close to the star). On 
the other hand, for cr < the sign is positive, which corre- 
sponds to heating. 

To calculate the function 7(i, cr, /i*), we use the numer- 
ical quadrature (KKID . First, the integral over x can be 
efficiently computed using a Hermite quadrature formula. 
Quadrature weights and knots were computed using the 
subroutine IQ PACK, which is an implem entation of a method 
described by iKautskv fc Elhavl (|l982f ) . A satisfactory ap- 
proximation can be obtained using 50 quadrature points. 
For large t (> 100), the Hermite quadrature formula be- 
comes inefficient and we use the simple trapezoidal rule. 
For angle integration, we used the Legendre quadrature 
formula with 5 quadrature points. Quadrature weights and 
knots were again compu ted using the subroutine IQPACK 
(IKautskv fc ElhavlfTgsa . 

To ensure the convergence of the model equations, we 
also included the derivatives of the GO heating term Eq. (|4]) 
with respect to the corresponding model variables, in the 
Newton-Raphson iteration step. 
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3. 

b 



0.03 



_ 0.02 



0.01 




Fig. 1. The dependence of 7(^,0-,^*) on t for selected spe- 
cific values of a and /z* 

3.2. The effect of GO heating 

To understand the role of GO heating in low density winds, 
we first neglect the frictional heating. This enables us to 
compare the magnitude of GO heating with the compet- 
ing cooling processes, namely radiative and adiabatic cool- 
ing. The cooling processes are included consistently in our 
NLTE models, although in the following we use optically 
thin radiative cooling after iRavmond et al.l ()1977t ). 



and adiabatic cooling 



J__d 

J.2 



(10) 



(11) 



where A(T) is the cooling function. For simplicity, we as- 
sume that riH ~ Uc, a? ~ 2fcT/mH. 



The importance of GO optically thick heating can be 
inferred by its comparison with other effects influencing 
temperature, namely adiabatic cooling and radiative heat- 
ing. Approximating the derivative in Eq. (fTTj) by using the 
velocity law v = i'oo(l — R*/^), assuming r ^ R^,, and 
vth.i/a ~ 1/4, then the ratio of GO heating and adiabatic 
cooling is from Eqs. ([TT|) 



GO 



ad 



47r 7th- 



E 



(13) 



where Voo is the wind terminal velocity, 7th is the value 
of j{t,a,fif) for t ^ 1 (typically according to Fig. [T] 
7th ~ 0.01), and M is the wind mass-loss rate. In scaled 
quantities, Eq. reads 



/ Vqo \ 
QP° .,no Vl03kmS-lJ 

gad -^■'^ ^ a \ 



IGOR, 



f ( M 

VlOkms-i/ 



E 



U015 



10-7 M© year- 1 



10~^ erg cm~^ 



(14) 



For stars with large radii (i?, > 10 Rq) and sufficiently 
high number of optically thick lines, GO heating may be 
comparable to the adiabatic cooling and may influence wind 
temperature. This was already found by GO, as Eq. ([Til) 
corresponds to Eq. (28) of GO. 

A comparison of GO heating given by Eq. (|4]) with ra- 
diative cooling given by Eq. (jlOp implies that 



GO 



Q 



rad 



1287r'^r'^u'^m^7th 
c2A(T)Af2 



thick 



(15) 



3.2.1. Low-density winds 

The size of the GO heating/cooling term may be inferred 
from the plot of the function 7(t, cr, /i*) in Fig. [T] For 
winds of low density, the line optical depths are small, 
(t/(l -I- (J^il) < 1), and we can use the Taylor expansion 
of the exponentials in Eqs. (O, ([7]). Assuming the line pro- 
file given by the Gaussian function, we derive for a > the 
approximation 



7(i,cr,^* 




arctan-v/CT — arctan (12) 



and in a similar way also for a < 0. Because it follows from 
Eqs. (HJ - (ini) that t ^ p, the GO heating/cooling term 
in Eq. ([T^ is proportional to p^. On the other hand, the 
adiabatic cooling term in Eq. (jllD is linearly proportional to 
density. Therefore, for very low wind densities when all lines 
become optically thin, adiabatic cooling begins to dominate 
and GO heating is not important. 

3.2.2. High-density winds 

From Fig. [1] it follows that for strong lines t ^ 1 the GO 
heating function j(t,a,fj,^) does not significantly depend 
on t. The value of (5p° for optically thick lines therefore 
basically depends on the number of these lines. 



or, as suming A(T) = 5 x 10-23(T/10^ K) (|Ravmond et al.l 
Il977f ) and approximating v Voo in scaled quantities 



GO 



")rad 





M 



T 



■ viooR©/ Vio^kms-iy Vio^'K 



-1/2 



10-7 Mq year- 



^ Vl0i5s-V 



10 erg cm ^ 



(16) 



which is similar in value to Eq. (31) of GO. If there is a 
sufhciently high number of optically thick lines, then GO 
heating may be comparable to the radiative heating in the 
outer regions of hot star winds. 

3.3. Numerical results 

An example of the numerical results obtained for the NLTE 
model M999-1 (see Tabled]) for different mass fractions of 
heavier elements Z is given in Fig. [2] The effect of the 
GO heating is the strongest for the model with the high- 
est metallicity (the highest mass- loss rate) , whereas for low 
metallicity models these effects are negligible. 

The trends displayed in Fig. [5] can be explained using 
the discussion presented in Sect. The dominant cooling 
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30000 




25000 



20000 



15000 



five-component model 
one-component model 



10000 



10 15 



20 



25 



Fig. 2. The effect of GO heating in the NLTE wind model 
M999-1 (see Table[T]) obtained for different mass fractions 
of heavier elements. Multicomponent effects were neglected 
here. 



mechanism for stars with a very low mass-loss rate (very 
low metallicity) is the adiabatic cooling (see Sect. I3.2.ip . 
For a star with Z = 10"^, the GO heating does not sig- 
nificantly influence the wind temperature. For stars with 
higher mass-loss rates (higher metallicity), the radiative 
and GO processes in addition to the adiabatic cooling in- 
fluence the wind temperature. For example for Z — 10^^ 
and r ~ lOi?*, we can derive froi n Eqs. (T Hl) and using 
corresponding wind parameters (jPaoer J ) that the ratio of 
GO heating to both adiabatic cooling and radiative heating 
is of the order of unity. The GO heating may therefore in- 
fluence wind temperature. For Z = 10"^, the mass- loss rate 
increases, but on the other hand both the terminal velocity 
and the number of optically thick lines increase, enhancing 
the importance of the GO heating. 

Although the GO heating may be important for temper- 
ature balance, in the calculations presented here it never 
heats the wind to temperatures significantly higher than 
the stellar effective temperature. The relative importance 
of this effect for temperature balance of the wind of stud- 
ied evolved stars is connected to their large radii and in 
some cases also with large wind terminal velocities. For 
present-day hot stars with smaller radii, lower wind ter- 
minal velocities, and iron as a significant wind driver the 
effect is less important, i.e., as iron lines become optically 
thin in outer wind regions, their contribu tion to GO heat- 
ing is likely to be relatively small (see also lVink et al.l[T999l 
iPuls et al.l[2000l) . 



4. Multicomponent models 

After discussing the role of GO heating separately, we 
present here detailed NLTE five-component wind models 
(described in Sect. 12. ip with all relevant heating/cooling ef- 
fects included, i.e., radiative, adiabatic, GO, and frictional 
ones. These models are supplemented by hydrodynamical 
simulations of two-component flow (see Sect. 



4.1. High density winds 

For stars with high density winds, the velocity differences 
between individual wind components are much smaller 
than corresponding mean thermal speed. For these stars. 





10000 




0.018 




0.016 
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0.006 




0.004 
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Fig. 3. Multicomponent effects in the NLTE wind model 
M999-1 (see Table[I]) for Z = lO'^. Top: Comparison 
of electron temperature in the five-component and one- 
component wind models. Bottom: Non-dimensional velocity 
difference (Eq. (jl8p ) between passive component p (hydro- 
gen and helium) and carbon, nitrogen, and oxygen in the 
five-component wind model. 



frictio nal heating is negligibl e and decoup ling doe s not 
occur (Castor et aD 19 76. Springmann fc Pauldrach [l992i 
iKrticka fc Kubatll200lL fVotruba et al.ll2007t) . The~winds of 
these stars can be adequately described by one-component 
models (e.g.. Fig. [3]). 



4.2. Winds with frictional heating 

For winds with lower densities, velocity differences between 
individual wind components become comparable to the 
mean thermal speed, and momentum transfer between the 
wind components becomes inefficie nt. This may cause fric- 
tional heating of the wind ( Spring mann fc Pauldrachlll992i 
ICurelll992l IKrticka fc K^atil2nnifl . 

In some cases, the wind temperature may increase to 
values of the order of 10^ K (see Fig. 2]) due to friction. 
Howeve r, the increase is typically smaller. Frictional heat- 
ing (see lBurgerslll969l) becomes important when it is com- 
parable to radiative cooling Eq. (ITU)) , i.e.. 



lnAG(a:pi) {v, - v^) = nnn,K{T), (17) 



where In A is the Coulomb logarithm, A(T) is the cooling 
function introduced in Eq. (jlOp . Tpj is the mean tempera- 
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Fig. 4. Frictional heating in the NLTE five-component 
wind model M200-1 (see Tabled]) for Z = IQ-^. Top: 
Temperature of individual wind components (oxygen, 
hydrogen and helium component, and free electrons). 
Temperature of carbon and nitrogen components is nearly 
the same as the hydrogen and helium one. Bottom: 
Nondimensional velocity difference between passive com- 
ponent p (hydrogen and helium) and nitrogen and oxygen. 



ture, and the dimensionless velocity difference 



2fc(miTp+mpTi) 
mimp 



(18) 



The subscripts denote values of the radial velocity v, num- 
ber density n, charge q, atomic mass to, and tempera- 
ture T of heavy ions (i) and passive component (p) . Using 
the approximations ric ~ tt-h, G(a;pi) « -^pi ~ 

Vi — Vp « axpi, approximating the passive component b; 
hydrogen, and using the formula in lRavmond et all ()197 
for the radiative cooling (see also Eq. ([TB]) '). we can ap 
proximately determine the temperature at which frictional 
heating is balanced to be 



I 



Tfric = 9600 K 



0.01 



2/3 



TOH 



2/3 



9i 



4/3 



ooT 



4/3 



We note that the a nondimensional velocity differ ence can 
be rou ghly calculated using Eqs. (18) or (23) of iKrtickal 
(|2006|) . If the wind temperature is much higher than Tfrjc, 
then frictional heating is negligible. On the other hand, 
if the wind temperature is equal or lower than Tfric, then 



Fig. 5. Calculated velocity (top) and density (bottom) pro- 
file from hydrodynamical simulations of two-component 
wind for the model MlOO-2. Heavier ions are denoted us- 
ing dashed line, hydrogen and helium using solid line. The 
components decouple at r 2 . 



frictional heating may infiuence wind temperature and in- 
crease it to the value of about Tfiic. Since many studied 
giant stars have winds that have xm of the order of 0.01, 
frictional heating may influence their wind temperature for 
the mass fraction of a given element Zi w 0.01. Moreover, 
the effect of frictional heating is insignificant for stars of 
very low metallicity Zi ^ 0.01. 



4.3. Decoupling of wind components 

For very low density winds, the frictional force becomes in- 
efficient e nabling the dynamical decoupli n g of wind com- 
pone i its (ISoringmann &: Pauldrach 1992L lOwocki fc Puis! 



l2002llKrticka fc KubatJl2002H Votruba et al.ll2007l^ ." In most 
cases, the decoupling is stimulated by frictional heating, 
which leads to an effective decrease in the frictional force 
(due to the dependence of frictional force on the tempera- 
ture). 

If decoupling occurs for velocities larger than the cor- 
responding escape speed, then the wind mass-loss rate 
remains basically unaffected by decoupling (see Fig. [S]). 
Here we plot results of hydrodynamical simulations of a 
two-component win d for the model MlOO-2 (see Tabled]). 
The C AK constants (Cast or. Abbott fc Kleinlll975llAbbottl 
Il982[ ). which characterize the effectivity of the momen- 
tum transfer from photons to absorbing ions, are a — 0.6, 
k = 0.027, and 6 = 0.1 (corresponding to Z = 5 x 10"^). At 
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Fig. 6. HR diagram for different types of stellar wind. 
Crosses denote studied stars, gray lines marked by eor- 
responding mass a re evolutionary tracks calculated by 
iMarigo et al.l ()2001[ ). solid lines denote approximate loca- 
tion of borders below which the fall back of hydrogen and 
helium occurs, and dashed lines denote approximate loca- 
tion of a border below which decoupling occurs (for a given 
mass fraction of heavier elements). 



some point, the absorbing ions begin to decouple from the 
passive plasma and accelerate rapidly (the so-called run- 
away effect). Because this occurs at velocities larger than 
the escape speed, all wind components finally leave the 
star. On the other hand, if decoupling were to occur at 
velocities smaller than the escape speed, then the hydrogen 
and helium components would be unab le to leave the star 
(jPorter fc Skouzalll999l iKrticka fc Kub at 2001). Fo r very 
low m etallicities, a purely metallic wind may exist (jBabell 
119951 luliglaub.2008lV 

However, the final fate of the decoupled material re- 
mains unclear, because the two-stream instabilities may 
change the nature of the solution. 

A HR diagram showing the different types of stellar 
wind is given in Fig. [S] For evolved stars of the highest lu- 
minosities, the multicomponent effects are important only 
for very low metallicities of the order oi Z k, 10~^. On the 
other hand, for less-massive stars decoupling can occur even 
for a mass fraction of heavier elements comparable to the 
solar one [Z « 10~^). The shape of border lines in Fig.|6]is 
given mainly by the mass-loss rate, and partly also by the 
terminal velocity, wind temperature, and charge. On aver- 
age, decoupling occurs for mass fraction of CNO elements 
lower than roughly 



For even lower metallicities, lower than 

L 



log^i 



back 



0.40- 0.65 log 

/ Tcs 
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VIO^K 



-3.1 



off 



105 K 



(20b) 



the passive component decouples from hydrogen and helium 
for velocities smaller than the escape speed and may fall 
back onto the stellar surface. 



5. Discussion 

5.1. Wind limits 

There exists a limiting mass-loss rate of heavier elements 
j^^max ]-,giow which hydrogen and helium remain in the stel- 
lar atmosphere and a purely metallic wind exists with a 
mass-loss rate Mj < M^^''''. As shown bv lUnglaubl (j2008[) . to 
achieve a hydrostatic hydrogen solution in the atmosphere, 
the magnitude of the gravitational acceleration should be 
larger than that of the hydrogen acceleration due to friction 
with heavier elements, i.e., 



GM nnni 47rg^gf 

> — lnAG(a;Hi) 



PH 



(21) 



Using the same approximations as in Eq. (llOp and the con- 
tinuity equation, we obtain the condition 



V 2 qllnk qf 



SxlO^^Mgyear 



.-1 



g2 InA 

TOi \ / M 

ma 



T 



IMr/IiO^K 



3/2 



3L 



-2 



(22) 



Since in most considered cases in this study the derived 
mass-loss rate of heavier elements is higher than M'"'''', 
hydrogen (and also helium, see Sect. I5.4p are driven out 
of the stellar atmosphere for the metallicities studied here. 
Subsequently, they may leave the star if there is no decou- 
pling or if decoupling occurs for velocities larger than the 
escape value. 

For metallicities lower than those studied here, the ra- 
diative force is, however, insufficiently strong to drive the 
wind containing the hydrogen and helium ions. At these 
very low metallicities, purely metallic winds may occur, if 
the radiative acceleration is large enough. The mass-loss 
rate of such a wind is lower than that given by Eq. ([2^ . 
Since this condition does not depend on density, we note 
that the value of the limiting metallic mass-loss rate M™'^^ 
is nearly constant throughout the stellar atmosphere (with 
some variations caused only by the temperature and charge 
variations) . 



log Zdoc 



0.48- 0.60 log 



L 

VLr. 



1 + 1.49 



105 K 



3.1 



105 K 



(20a) 



5.2. Mass-loss rate and evolutionary calculations 

From th is and previous studies ([Krticka fc KubatI I2006L 
iPaper J ) , we obtain the following picture of the stellar winds 
of hot first stars with an initially pure hydrogen-helium 
composition. 
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Pure hydrogen-helium stars do not have any line-driven 
wind. Stars very close to the Eddington limit with F > 
0.859 can have very weak (M < 10^^** year" ^) pure 
hydrogen winds due to light scattering on free electrons 
(jKrticka fc Kubatll2006[ ). As a result, the influence of winds 
on the evolution of these stars is negligible. In rotating 
stars, mass loss by means of a decretion disk (Lee et 
al. 1991, Ekstrom et al. 2008) may be of some impor- 
tance. Furthermore, massive stars may lo se mass by means 
of either a ?7 Car type of explosions (|Vink fc de Koteil 
L2005, Smith & Owocki 2006) or a super-Eddington outflow 
(jOwocki et al.i 120041) . This outflow may more easily exist 
in evolved stars as the nuclearly processed core has a lower 
number of free electrons per nucleon than the hydrogen- rich 
envelope. 

As soon as heavier elements are synthesised in the stellar 
core and transported t o the stellar surface, a pure metallic 
wind may be produced ()Babellll9"95l ) if the radiative force on 
metals is large enough. These very weak stellar winds have 
mass-loss rates lower than those given by Eq. (|22|) (see also 
lUngl aub 2008) . Although this type of outflow probably does 
not influence stellar evolution, it may in fluence the strat- 
ificat ion of the stellar atmosphere (e.g., iLandstreet et al.l 
[19981) . 

For metallicities higher than those given by Eq. (18) of 
iPaoer J . wind that also contains hydrogen and helium may 
exist. For metallicities lower than those given by Eq. (I20bp . 
the wind is very weak, wind decoupling occurs at velocities 
lower than the escape value, and hydrogen and helium may 
fall back on the stellar surface. The final fate of such a wind 
is however unclear because realistic simulations describing 
their behavior are not yet available. 

For metallicities higher than those given by Eq. (I20bp . 
either winds are denser and decoupling occurs at velocities 
larger than the escape value or there is no decoupling at all. 
The wind mass-los s rate is g iven by the formulae of, e.g., 
IVink et al.l (|200lD . iKudritzkj (|2002D . or lPaoer 1 

5.3. High velocity particles in primordial minihaloes 

After decoupling, the CNO particles may be accelerated to 
velocities of the order of u w lO'^kms"^, i.e., about 0.1c. 
The energies of these particles are of the order of 1 MeV, 
and the consequent decoupling of the wind components pro- 
duces the first low-energy cosmic-ray particles. The slowing- 
down time of these particles in the primordial minihaloes 
can be considerably large. If the interstellar material were 
alread y ionised, th e slowing-down time would be estimated 
to be ((Diveijl200H) 



mjTOH 



Anqfq^nii In A 



(23) 



Assuming that the typical hydrogen densi ties inside mini- 
haloe s are of the order of lO^cm"'^ (e.g., iMachacek et al.l 
l200lh . the slowing-down time is of the order of 10^ years. 
Consequently, rapidly moving CNO particles may easily 
penetrate the halo and influence its chemical composition. 

The aforementioned analysis is relevant if the influence 
of fast metallic wind on the interstellar medium can be 
neglected. In the opposite case, if acceleration of interstel- 
lar medium is non- negligible, the medium is swept up by 
a metallic wind, and its density increases. As a result, the 
slowing down time (j23p shortens, the coupling between met- 
als and interstellar medium strengthens, and the process 



repeats itself until a structure typica l of wind-driven bub- 
bles develops fe.g.. iFrever et aLll2003[ ). To assess the likeli- 
hood of this scenario, we calculate the ionized interstellar- 
hydrogen accel eration caused by its friction with the metal- 
lic wind to be (|Burgerslll969() 



Airqfq^ni 

gmc = r^^\nAG{xm) 

mukTm 



gfggtAfjlnA 

2 9 ' 

rriim^r'^V'^ 



(24) 



where we used the approximate form of G'(a;Hi) ~ l/(2a;Hi) 
valid for xm ^ 1 and the continuity equation of metals. 
Writing the metallic mass-loss rate as Mi — (M'^^^ (see 
Eq. ) , the frictional acceleration is 



(25) 



where we have neglected the difference in charges between 
atmosphere and interstellar medium, and g(r) is the magni- 
tude of the gravity acceleration at a given point. Although 
may be of the order of 10^, the ratio of the thermal speed in 
the atmosphere to the metallic wind velocity is a/w « 10"'^. 
Hence, the frictional acceleration of the interstellar medium 
is seven orders of magnitude lower than the gravitational 
acceleration and does not signiflcantly influence the dynam- 
ics of the interstellar medium. 

For homogeneous winds (i.e., those also including hydro- 
gen), both the mass-loss rate and the ratio a/v are larger. 
In this case, the frictional force influences the dynamics of 
the interstellar medium, the slowing-down time of Eq. (1231) 
is short, and the problem becomes a hydrodynamical one. 

5.4. He-free winds of cool stars? 

Being heavier and less fragile than neutral hydrogen, neu- 
tral helium may remain in the stellar atmosphere, while 
hydrogen flows with heavier elements in the stellar wind. 
Taking into account that collisions with hydrogen are more 
important for helium acceleration than collisions with heav- 
ier ele ments, helium re mains in the stellar atmosphere if 
(Krticka fc Kubatll2004 c.f. Eq. ^) 



GM 



9„ /cTr 



> 



In AG(a;HQ), 



(26) 



where the subscript a denotes quantities corresponding to 
helium. Using the approximate formula for G(xHa) and the 
hydrogen continuity equation, we obtain a condition for the 
hydrogen mass-loss rate 



Mh < 3 



g2 InA 



10"^®Moyear~i 



M 



( ^ 



VlO^K 



3/2 



(27) 



From this, it seems that if the helium charge qa is very low, 
i.e., if most of the helium atoms are neutral, then helium 
may remain in the stellar atmosphere, while hydrogen and 
heavier elements produce a wind. The other possibility of 
helium- free wind, i.e., a very low hydrogen mass- loss rate 
-^H ^ 10~^^ M0 year~^, is not possible here because all 
studied stars have mass-loss rates higher than this value. 

However, if helium atoms are neutral, collisions between 
neutral helium atoms and ionised hydrogen can accelerate 
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helium into the stellar wind. Neutral helium remains in 
the stellar atmosphere if the magnitude of the frictional 
acceleration is smaller than the magnitude of the gravity 
force, i.e., 



(28) 



where fina — ^nY^nia/ {ttiy^ + ma) is the reduced 
mass, and (cru)Ha is the momentum transfer rate co- 
efRci ent for coll isions betwe en neutra l helium and pro- 
tons (jPinto fc Ga Ui 2008, Krs tic fc Shu hz 1999). Adopting 
(tTw)Ha ~ lO^^cm'^s^^ and using the continuity equation, 
the condition given by Eq. (j28p can be rewritten as 



2 X 10"^^ Mq year" 



M 



IM 



-0 



(29) 



For cooler massive stars, collisions between neutral helium 
and protons alone are able to accelerate helium if the wind 
is strong enough. 

In models of cooler stars that we have studied, there is 
indeed a large fraction of neutral helium, but the charge of 
ionized helium is always high enough to accelerate helium 
into the wind (see Eq. (I?71) '). For stars cooler than those 
studied here (with T^s < 10** K), collisions of neutral helium 
with protons dominate the acceleration of helium in the 
stellar wind. 

5.5. Charging of first stars and first magnetic fields 

The problem of the escape of electrons from the s tellar at- 
mosp here was considered by Milne (192^, see also I Johnsod 
Il925f) . He demonstrated that electrons, due to the radiative 
forces acting on them, can escape from the star. However, 
he concluded that the resulting positive charge of the star 
would soon prevent any additional loss. The problem of 
charging the first stars by means of their winds remains 
an interesting one because of its possible implications for 
generation of the first magnetic fields. 

Our models of the multicomponent stellar wind do not 
allow us to check the possibility of a star increasing its 
charge, since the zero current condition was used to deter- 
mine the electron velocity at the wind base. The possible 
inclusion of an electron regularity condition corrupts (al- 
though only very slightly) the zero current condition at the 
upper boundary and can lead to the charging of the star. 

However, the physical significance of both the electron 
critical point and the regularity conditions is questionable. 
From the hydrodynamic point of view, these conditions can 
only be used if they correspond to the point where the prop- 
agation speed of some type of disturbance is met. In gen- 
eral, the characteristic equations for the multicomponent 
flow may n ot necessarily describe the dissemination of the 
waves fe.g.- lKrticka fc Kubatll200l . Thus, hydrodynamical 
simulations (or at least a linear analysis of hydrodynamical 
equations) should be used to prove the significance of the 
electron critical point condition. Only an analysis of this 
type will be able to answer the question of whether stellar 
charging via stellar winds is conceivable. However, if the 
possibility of stellar charging via stellar winds were proven, 
it could lead to the creation of first magnetic fields in the 
Universe. 



5.6. The effect of wind inhomogeneities 

We have modelled hot star winds by neglecting small- 
scale inhomogeneities . These in homogeneities are expected 
on both theoretical (jOwocki e t al. 1988, Feldmeier et all 
'199 7D and observational grounds (jBouret et al.l l2003l 
iMartins et al.ll2005l iPuls et al.l[2006D . However, we expect 
that even in structured winds the inefficient transfer of mo- 
mentum between the wind components may strongly affect 
the wind structure. We provide at least an estimate of the 
wind parameters for which these problems may occur and 
we outline the inclusion of these effects in the evolutionary 
calculations. 



6. Conclusions 

We have studied the effect of the Gayley-Owocki (Doppler) 
heating and multicomponent flow structure in CNO driven 
winds of hot stars. The parameters of these stars were se- 
lected to represent massive initially pure hydrogen-helium 
(Pop HI) stars. 

For the flrst time, we have included the GO heating 
term directly using atomic linelist and NLTE calculations. 
We have shown that GO heating is important especially for 
winds of CNO enriched flrst stars with high metallicities 
{Z « 0.01). In these winds, the GO heating can compete 
with radiative cooling because both the number of strong 
lines and the wind terminal velocity are large. On the other 
hand, for stars with low metallicities {Z < 0.001) there 
are an insufficient number of strong lines and the adiabatic 
cooling dominates. 

The effects of multicomponent flows are important espe- 
cially at low metallicities {Z < 0.001) in the case of evolved 
stars, and for relatively high metallicities (Z w 0.01) for 
main-sequence stars. The frictional heating itself does not 
influence the wind mass-loss rate. On the other hand, de- 
coupling probably leads to a zero mass- loss rate of hydrogen 
and helium if it occurs at velocities lower than the escape 
one. We have developed an approximate formula that esti- 
mates the minimum metallicity above which hydrogen and 
helium leave the star. 

The decoupling of radiatively accelerated metals from 
hydrogen and helium leads to generation of particles with 
typical energies of the order of IMeV, i.e., the first stars 
may be the first sources of low-energy cosmic rays. We have 
shown that these particles easily penetrate the interstellar 
medium of a given minihalo. We have discussed the pos- 
sibility of charging of first stars via their multicomponent 
winds. 

Wind models presented here can also be used to describe 
the winds of possible subsequent generations of CNO rich 
stars and low-luminosity stars of solar chemical composi- 
tion. 
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